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TRANSIENT  HYDROTHERMAL  STRESSES  IN  COMPOSITES: 

COUPLING  OF  MOISTURE  AND  HEAT  WITH  TEMPERATURE  VARYING  OIFFUSIVITY 

by 

G.  C.  Sih  and  M.  T.  Shih 
Institute  of  Fracture  and  Solid  Mechanics 
Lehigh  University 
Bethlehem,  Pennsylvania  18015 

ABSTRACT 

In  this  paper,  the  influence  of  coupled  diffusion  of  heat  and  moisture  on 
the  transient  stresses  in  a composite  is  investigated  analytically  where  the 
moisture  diffusion  coefficient  is  taken  to  be  temperature  dependent  while  the 
thermal  diffusion  coefficient  is  kept  constant.  There  is  no  a priori  reason 
why  moisture  and  temperature  should  be  uncoupled  such  that  each  will  obey  the 
simple  diffusion  theory,  particularly  without  reference  made  to  the  initial  and 
boundary  conditions  of  a particular  situation.  A study  of  the  coupled  diffusion 
equations  were  made  by  a finite-difference  scheme  allowing  for  time-dependent 
changes  in  the  humidity  and  temperature  of  the  environment.  The  appropriate 
transient  boundary  conditions  are  specified  on  the  surfaces  of  an  infinite  plate 
Numerical  calculations  were  carried  out  for  the  T300/5208  graphite  fiber-rein- 
forced epoxy  matrix  composite  in  which  the  nonuniformity  of  moisture  and  temper- 
ature is  evaluated  for  sudden  changes  in  the  surface  moisture  and/or  temperature 
The  coupling  effect  between  temperature  and  moisture  is  found  to  be  most  signifi 
cant  when  the  plate  undergoes  a sudden  change  in  surface  temperature  while  the 
surface  moisture  concentration  is  held  constant.  The  present  findings  indicate 
that  the  stresses  due  to  coupling  can  deviate  from  the  uncoupled  results  any- 
where from  20  to  80  percent  depending  on  the  surface  temperature  gradient.  This 
suggests  the  need  to  perform  additional  experiments  for  evaluating  the  coupled 


diffusion  phenomenon  and  its  influence  on  the  mechanical  behavior  of  epoxy-resin- 
composites. 


i 


INTRODUCTION 

Absorption  of  moisture  by  composites  causes  dimensional  changes  through  non- 
uniform  expansion  and/or  contraction  of  material  elements  which  in  turn  leads 
to  internal  stresses  and/or  strains.  If  the  physical  process  is  of  a more  ac- 
tive type,  caused  by  capillary  flow  of  fluid  into  voids,  stresses  can  occur  even 
when  the  macroscopic  deformation  is  uniform.  These  fluid- induced  stresses  could 
explain  the  lowered  stiffness  and  strength  of  composites.  Moreover,  the  thermal 
environment  may  also  interact  with  moisture.  For  example,  voids  and  microscopic 
cracks  open  as  temperature  is  increased,  and  more  fluid  is  absorbed  into  the 
openings;  the  trapped  fluid  causes  growth  of  the  flaws  when  temperature  is  re- 
duced suddenly.  Subsequently,  the  material  is  capable  of  absorbing  more  moisture 
more  quickly  than  before  the  thermal  gradient  were  applied.  Such  a process,  if 
continued,  could  lead  to  eventual  failure  of  the  composite. 


There  is  an  urgent  need  to  quantitatively  assess  the  moisture  and  temperature 
interaction  effects  which  can  affect  the  physical  properties  [1]  and  stress  char- 
acteristics [2,3]  of  epoxy- composite  materials.  At  present,  the  composite  tech- 
nology community  is  very  much  concerned  with  this  problem  because  more  and  more 
of  the  composites  with  epoxy-resin  matrix  are  being  used  for  constructing  air- 
craft structures.  No  confidence  level  in  design  could  be  established  unless  the 
behavior  of  these  materials  in  the  presence  of  adverse  environments*  is  under- 
stood. 


? 

The  influence  of  chemical  reaction,  which  may  also  play  an  Important  role  in 
many  instances,  will  not  be  addressed  in  this  communication. 


. 

As  mentioned  earlier,  the  diffusion  of  moisture  and  temperature  is  intimately 
related  to  the  degradation  of  both  the  strength  and  stiffness  of  composites. 

Increasing  the  temperature  and  moisture  content  in  the  composite  can  generally 
accelerate  the  degradation  process.  Previous  investigations  [1,2]  have  assumed 
as  a priori  that  these  two  effects  are  uncoupled.  Such  a decision,  however, 
cannot  be  made  intuitively  but  should  rely  on  an  analysis  of  the  coupling  theory 
for  the  particular  problem  involving  specific  initial  and  boundary  conditions. 
Information  on  this  subject  is  to  say  the  least  very  lacking  at  this  time  and  no 
general  conclusion  can  be  given.  The  significant  variables  in  such  a study  must 
at  least  involve  time,  relative  humidity  of  the  environment,  temperature,  rele- 
vant physical  constants,  etc.  The  realistic  prediction  of  the  moisture  and  tem- 
perature distribution  in  composites  depends,  of  course,  on  selecting  the  appro- 
priate governing  equations.  To  this  end,  various  phenomenological  arguments 
were  presented  in  [4]  to  explore  a variety  of  coupled  equations  describing  the 
simultaneous  diffusion  of  moisture  and  heat.  Five  different  physical  models 
were  obtained.  The  coefficients  in  these  models  are  associated  with  the  basic 
thermodynamic  properties  of  the  solid  and  can  be  related  to  one  another.  This 
was  made  possible  because  the  basic  form  of  the  coupled  equations  for  the  five 
models  turned  out  to  be  the  same. 

The  purpose  of  this  investigation  is  to  develop  an  analytical  model  and 
technique  for  calculating  the  nonuniform  moisture,  temperature  and  stresses  in 
composite  systems.  The  T300/5208  graphite/epoxy  system  was  used  in  the  numeri- 
cal calculation,  because  diffusion  data  as  well  as  the  variation  of  the  moisture 
diffusion  coefficient  with  temperature  for  this  material  are  available  [5].  A 
finite  difference  computer  program  was  developed  for  solving  the  coupled  diffu- 
sion equations  with  transient  boundary  conditions  on  moisture  and/or  temperature. 
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Examples  and  numerical  calculations  are  provided  for  moisture  and/or  temperature 
diffusing  into  a plate  from  its  surfaces.  The  plate  is  initially  at  a uniform 
temperature  with  a given  moisture  content  distributed  uniformly  throughout  the 
plate.  Suddenly,  the  temperature  and/or  moisture  at  the  plate  surfaces  are 
changed  and  maintained  constant  thereafter.  The  corresponding  stresses  are  also 
calculated  as  a function  of  time  while  the  numerical  results  for  other  quanti- 
ties of  interest  are  displayed  graphically. 

COUPLED  DIFFUSION  EQUATIONS 

The  derivation  of  the  basic  equations  coupling  heat  and  moisture  content  in 
a solid  has  already  been  discussed  in  [4]  and  will  not  be  repeated  here.  The 
second  model  in  [4]  refers  to  the  equations 


DV2C  - (C-AT)  = 0 

m2T  - (T-vC)  = 0 


0) 


where  V2  is  the  Laplacian  operator  in  the  space  variables.  In  equations  (1),  T 
is  temperature  and  C is  the  mass  of  moisture  per  unit  volume  of  void  space  in 
the  solid.  The  diffusion  coefficients  D and  V have  units  of  area  per  unit  time, 
and  the  coupling  coefficients  A and  v have  units  of  mass  per  unit  volume  per 
unit  temperature  and  the  reciprocal,  respectively.  These  equations  are  rela- 
tively easy  to  solve  when  the  coefficients  are  constant  and  boundary  values  of 
temperature  and  moisture  content  are  held  constant  between  occasional  moments 
of  sudden  changes  [3]. 

Solutions  of  equations  (1)  can  be  examined  to  determine  the  influence  of 

coupling  on  the  overall  behavior  of  composites.  It  has  been  noted  experimentally 
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that  the  moisture  diffusion  coefficient,  0,  depends  on  temperature  by  a relation 
of  the  form  [6] 


D « Dq  exp(.-E0/RT) 


U) 


in  which  EQ  is  the  energy  required  for  one  unit  of  mass  to  move  into  the  solid, 

R is  the  gas  constant,  and  T is  the  absolute  temperature.  The  diffusion  coeffi- 
cient, D,  is  usually  measured  in  an  experiment  such  as  that  described  in  [1]  for 
the  uncoupled  diffusion  equation.  The  meaning  of  D is  altered  when  coupling  is 
present  and  the  values  of  D in  [1]  can  differ  by  an  amount  which  depends  on  the 
values  of  V and  Xv  [7].  When  D is  temperature  dependent  as  shown  in  equation  (2), 
then  the  following  coupled  theory  must  be  used: 


V • (DVC)  - (C-XT)  = 0 
W2T  - 1^-  (T-vC)  = 0 


(3) 


where  V is  constant  throughout  the  present  discussion.  When  D is  a function  of 
temperature,  equations  (3)  are  nonlinear  and  a numerical  scheme  for  solving  equa 
tions  (3)  is  required. 

For  the  problem  at  hand,  only  the  moisture  and  temperature  changes  in  the 
plate  thickness  or  z-direction,  Figure  1,  is  considered  and  hence  V2  = 3i/3z2. 

It  is  expedient  to  introduce  the  dimensionless  space  and  time  variables 


(4) 
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in  which  h stands  for  the  plate  thickness.  Equations  (3)  can  thus  be  expressed 
in  terms  of  £ and  0 as 


In  what  follows,  two  types  of  transient  boundary  conditions  will  be  treated, 
namely  sudden  change  in  moisture  and  temperature. 

SUDDEN  CHANGE  IN  MOISTURE 

Consider  the  problem  of  diffusion  into  an  infinite  plate  as  shown  in  Figure 
1.  The  temperature  and  moisture  concentration  are  initially  uniform  at  the 
values  T..  and  C.,  respectively.  At  time  t=0,  the  moisture  at  both  surfaces 
z = ±h/2  are  suddenly  changed  to  C^.,  and  maintained  constant  thereafter.  The 
surface  temperature  of  the  plate  is  always  kept  at  T.. . These  conditions  may  be 
stated  as 

T(z,o)  = T.,  C(z,o)  = Ci  (6) 

and 

T(±h/2,t)  = T.,  C(±h/2,t)  = Cf  (7) 

for  time  t>0.  In  terms  of  the  nondimensional  variables  C and  0 in  equations 
(4),  the  solution  for  moisture  and  temperature  may  be  expressed  in  the  forms 
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(8) 


c(t)  = C1  + (Cf-C.)  f(£,8) 


T(t)  = Ti  + v(Cf-C.)  g(C,9) 


in  which  f(5,9)  and  g(5,9)  are  functions  t be  determined  from  the  conditions 
in  equations  (6)  and  (7).  Substituting  equations  (8)  into  (5)  yields 


3f  _ 1 rXv  3fg  . cra2f  , AB  3f  3q-,, 

‘ TTYTT  T?r  + ( ] +Ra ) * 3£-*^ 


36  ^-\x>  luQ  3^  L35 


35  35J 


3q  _ 1 rl  S^g  cr32f  j AB  3f  3q-ti 

#'T^{ir#+  F[35t  [T+lgP  35  35J} 


where  u„  stands  for 
o 


uo  = IT 


and  F is  given  by 


F * exp[-  T+BglT^T1 


The  quantities  A and  B are  defined  as 


Eo  . ',<cf-ci> 


A = RT7’  B ■ 


finite.  di^eAZnce.  mzthod.  Since  equations  (9)  cannot  be  solved  analytically, 
it  is  necessary  to  resort  to  approximate  numerical  methods.  The  method  of  finite 
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difference  is  adopted  to  replace  the  governing  partial  differential  equations 
and  the  associated  transient  boundary  conditions  by  the  corresponding  finite- 
difference  equations.  This  then  reduces  the  problem  to  a set  of  simultaneous 
algebraic  equations  which  can  be  easily  solved.  Referring  to  the  space  and 
time  interval  in  Figure  2,  the  first  of  equations  (9)  may  be  written  in  differ- 
ence form  as 


^m,n+l“^m,n  _ 1 rX\>  /^m+1  ,n'^^m,n+9m-l  ,nN  , r, \ rFm+l  ,n~^m,n+fm-l  ,n 

A0 { ^ ( [ATP ' (m,n)  [ [AC? 

, AB  /m+l  ,n~fm,n  w9m+l  ,n~9m,n^i 


while  the  second  of  equations  (9)  becomes 


3m,n+l"9m,n  _ 1 rl  ,9m+l ,n”29m,n+9m-l ,n,  . c/  , r*m+l  ,n’^m,n+^m-l  ,n 

A9 NXv  {u  ( JMV } F(m,n)  C ftcF 

0 

AB  /^tH-I  ,n"fm,nw9m+l  ,n~9m,nN-i-, 

n+Bg;jny7  ( — ac — )( — ac — )]}  (1< 


In  order  to  achieve  acceptable  accuracy  in  the  finite  difference  calculations, 
the  grid  size  in  space,  Az,  and  time,  At,  must  be  sufficiently  small  and  satisfy 
the  stability  requirement  that 


At  < exp(E0/RT) 


The  boundary  conditions  in  equations  (6)  and  (7)  may  then  be  written  in  terms 
of  f(C>9)  and  g(£,9).  They  become 


f(C.o)  = 0,  g(C,o)  = 0 
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f(±l,0)  = 1,  g(±l,9)  = 0 


for  9>0. 


A computer  program  was  developed  to  solve  equations  (13)  through  (17)  for 
the  functions  f(£,9)  and  gU,8)  from  which  the  moisture  and  temperature  through- 
out the  solid  can  be  determined. 

Ave/iage  mclstiLtz  quantLtieA.  In  this  model  [4],  the  mass  of  moisture  con- 
tained in  the  solid  per  unit  mass  of  solid,  M,  is  assumed  to  be  linearly  related 
to  the  mass  of  moisture  contained  in  the  voids  per  unit  volume  of  void  space,  C, 
and  the  temperature,  i.e.. 


M = aC  - ojT  + constant 


where  a and  oj  are  constants.  The  mass  of  moisture  contained  in  the  volume  of 
the  composite  per  unit  mass  of  solid,  m,  is  given  by 


m = - C + M = oj[y  - T]  + constant 

P A 


in  which  v is  the  fraction  of  unit  volume  of  the  voids  and  p is  (l-v)p  where 
Ps  is  the  density  of  the  solid  when  no  voids  are  present.  The  average  values 
of  these  moisture  content  quantities  are  defined  as 


C * ~ / CdV,  ft  * l / MdV,  m = l j mdV 
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The  total  moisture  in  the  voids,  solid  and  composite  are,  respectively,  vVC, 
pVfl  and  pVm. 


Now,  let  the  average  values  of  T,  C and  m be  defined  by  the  integrals 


, h/2 

T(t)  = Tu  / T(z,t)dz 

" -h/2 

, h/2 

C(t)=J  / C(z,t)dz  (21) 

n -h/2 

, h/2 

m(  t)  = r-  / m(z,t)dz 

" -h/2 


In  terms  of  f(£,0)  and  g(£,6),  equations  (21)  become 


1 

T(t)  - Ti  = J (Cf-C.)  / gU,9)dC 
, 1 

C(t)  - Ci  = \ (Cf-C.)  / f(5.9)d5  (22) 

1 

m(t)  - m.  = ^ (Cf-C.)  j [fU,9)  - AvgU,0)]dS 


In  view  of  equations  (6)  and  (7),  the  third  of  equations  (22)  may  be  put  into 
the  dimensionless  form  (m(t)-m.j)/(m^-m. ) which,  when  approximated  by  Simpson's 
rule  for  a fixed  time  0Q,  gives 


■ l (f>  <[f(l,90)  * 4f(2.0o)  + 2f(3,90)  * ...  - 4f(n-l,60) 


+ f(n,0Q)]  - Av[g(l,0o)  + 4g(2,9Q)  + ...  - 4g(n-l,eQ) 


+ g(n,0o)]}  (23) 
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Having  completed  the  preliminaries,  a numerical  example  will  now  follow. 


HumeAxcaZ  example..  Numerical  calculations  are  made  for  a T300/5208  epoxy 

2 

resin  plate  with  thickness  h = 0.2  cm.  The  constants  DQ  = 1.53  x 10  cm2/hr  and 
4 

Eq  = 1.25  x 10  cal/g-mole  are  obtained  from  [5J.  For  the  coupled  diffusion 

★ 

problem,  the  particular  values  of  u = 0.1,  A = 0.5  and  v = 0.5  are  chosen  for 
the  calculation  while  w is  an  arbitrary  constant.  The  choice  of  this  selection 
will  become  evident  subsequently.  Note  that  u = D/0,  should  be  distinguished 
from  uQ  in  equation  (10).  The  constants  A and  B in  equation  (12)  is  determined 
from  an  initial  temperature  of  T^  = 21 °C  = 294°K,  Ci  = 0 and  a gas  constant  of 
R = 1.986.  Hence 


A = 


o 

*T7 


vC. 


1.25  x 10* 
1.986  (294) 


= 21.4 


(24) 


-^=  1-7  x 10'3  Cf 


in  which  C^,  the  equilibrium  moisture  concentration,  can  be  obtained  from 


Cr  = 


OV' 


f • 


(25) 


for  different  relative  humidity  or  RH  of  the  environment.  In  equation  (25),  ys 
is  the  specific  humidity  measured  in  grams  of  water  per  lb  of  dry  air  and  p is 
the  density  of  the  ambient  air  in  units  of  g/cm3.  Another  important  quantity  in 
the  diffusion  analysis  is  the  relationship  between  equilibrium  moisture  content 


The  choice  of  this  selection  is  based  on  comparing  the  analytical  prediction  of 
the  percent  moisture  content  as  a function  of  /t  with  the  experimental  data  in 
[5],  Figure  3. 
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m^)  of  the  composite  and  RH  of  the  environment.  For  the  post-cured  T300/5208 
epoxy  resin,  the  relation  [5] 

m(»)  = 0.0155  (RH)  (26) 

may  be  used  in  which  RH  is  expressed  in  percent.  There  remains  the  appropriate 
selection  of  the  time  and  space  interval  before  carrying  out  the  finite  differ- 
ence calculations.  The  plate  in  the  z-direction  is  divided  into  seven  segments 
and  hence  Az  = h/7  while  At  must  satisfy  the  stability  condition  in  equation 
(15).  Results  are  expressed  in  terms  of  percent  moisture  content  m(t)  as  mani- 
fested by  the  weight  gain  of  the  composite: 

w ( t )-w . 

m(t)  = — ' — x 100  (27) 

Wi 

where  w(t)  is  weight  of  the  specimen  at  time  t and  w^  the  initial  dry  weight  of 
the  specimen. 

Figure  3 gives  a plot  of  m(t)  versus  /t  for  different  relative  humidities  of 
RH  = 13,  33,  52,  75  and  100".  The  moisture  diffusion  coefficient  D is  assumed 
to  be  temperature  dependent.  The  dotted  curves  represent  solutions  for  the  un- 
coupled theory  in  which  A = & = S = 0,  \v  = 0 and  u can  be  arbitrary.  They  dif- 
fer very  little  from  the  curves  for  the  coupled  theory.  The  values  of  u = 0.1, 

X = 0.5  and  v = 0.5  are  selected  such  that  the  coupled  solutions  in  Figure  3 
give  the  best  fit  to  the  experimental  data  in  [5].  Similar  results  can  also  be 
obtained  for  T^  = 43°C,  63°C  and  82°C.  Figures  4 to  8 show  the  variations  of 
moisture  content  m(t)  with  the  normalized  thickness  coordinate  2z/h  for  RH  = 13, 
33,  52,  75  and  1 00fo . Initially,  i.e.,  for  small  time  t,  only  the  region  close 
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to  the  plate  surface  experiences  moisture  while  the  center  region  of  the  plate 
is  not  affected.  As  time  increases,  moisture  is  penetrated  into  all  the  mate- 
rial elements  with  the  minimum  influence  at  2=0.  The  difference  of  m(t)  between 
z=0  and  z = ±h/2  increases  with  increasing  RH.  The  effect  of  initial  tempera- 
ture on  the  penetration  of  moisture  is  shown  in  Figure  9 for  a sudden  change  of 
RH  from  0"  to  100%.  Coupling  is  neglected  and  D is  taken  to  be  a constant. 

The  time  at  which  the  plate  reaches  moisture  equilibrium  is  seen  to  decrease  as 
T.  is  increased  for  a fixed  value  of  h = 0.2  cm. 

In  what  follows,  the  effect  of  sudden  change  of  surface  temperature  on  the 
moisture  uptake  in  the  plate  will  be  studied. 

SUDDEN  CHANGE  IN  TEMPERATURE 

In  order  to  study  the  influence  of  a sudden  temperature  on  a composite  that 
is  kept  in  a constant  moisture  environment,  it  is  necessary  to  use  the  coupled 
theory  since  the  moisture  diffusion  equation  for  the  uncoupled  theory  applies  . 
only  to  the  case  of  constant  temperature. 

Suppose  that  the  surface  temperature  on  the  plate  in  Figure  1 is  raised  from 
an  initial  value  of  T.  to  a final  value  T^  and  the  moisture  concentrations  at 
z = ±h/2  are  kept  constant  at  all  time.  Then,  in  addition  to  equations  (6),  the 
following  conditions  must  also  prevail: 

T(±h/2,t)  = Tf , C(±h/2,t)  = C.  (28) 

The  form  of  the  solution  expressed  in  terms  of  the  variables  t.  and  0 defined  in 
equations  (4)  is 
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(29) 


C(t)  = C.  + A(Tf-T.)  f(5,0) 

T(t)  - Ti  + (Tf-T.)  g(C,9) 

Without  going  into  details,  the  governing  differential  equations  become 


3f 

39 


32f 


AB 

rFgP 


if.  , l a2q 
3C  35J  UQ  W 


(30) 


la 

39 


1+A.v  32q 

u0  ^ 


AB  3f  la-, 
(1+Bg)2  3£  3?J 


where  uQ  is  given  by  equation  (10).  Refer  to  equations  (11)  and  (12)  for  the 

★ 

definition  of  A,  B and  the  function  F.  Making  use  of  equations  (29),  the  condi- 
tions in  equations  (2B)  may  be  written  as 

f(±l,9)  = 0,  g (±1 , 9 ) = 1 for  9>0  (31) 

As  in  the  previous  example,  equations  (30)  will  be  solved  numerically  by  the 
finite  difference  method. 

Finite  di^exenee  equation 6.  Equations  t30)  will  now  be  cast  into  the  fi- 
nite difference  form.  With  the  nondimensional  time  and  space  interval  as  chosen 
in  Figure  2,  the  following  expressions  are  obtained: 

fm,n+l‘fm,n  _ F(m,n)  r^l  ,n‘^m,n+^m-l  ,n  . AB  ,^m+l  ,n~fm,n  x 

as  • -ter [ rstr  o*69Bin)' ( — * 

« folA")]  + J-  (32) 

★ 

With  the  exception  of  B which  in  the  case  of  sudden  temperature  change  should 
read  as  (T^-TJ/T,  instead  of  that  in  equation  (12)  for  the  case  of  sudden  mois- 
ture change. 
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The  stability  requirement  for  selecting  the  relative  size  of  At  and  Az  in  the 


numerical  calculation  is  the  same  as  that  stated  in  equation  (15). 

Mo-UtuA&  content.  Following  the  definitions  of  the  various  moisture  param- 
eters as  discussed  earlier,  the  average  moisture  content  in  the  composite  per 
uni t mass  of  sol  id  is 

1 

m(t)  - m.  = ^ (Tf-Ti ) / [f(5,9)  - gU.9)]d£  (34) 

From  equation  (19),  it  can  be  shown  that 


mf  - mi  = - (jj(T^-T^  ) (35) 

which  when  substituted  into  equation  (34)  yields 

_ 

m(t)-m,  . , . 

1^L  - \ if)  ([gO,90)  + 4g(2,0o)  ♦ ...  ♦ g(n,90)] 

- Cf(1.0o)  + 4f(2,0o)  + ...  + f(n,e0)]}  (36) 
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Simpson's  rule  has  been  applied  for  evaluating  the  integral  in  equation  (34)  at 

0 - e0. 

Numerical  exempli..  Referring  to  the  conditions  prescribed  by  equations  (6) 
and  (28),  the  moisture  concentrations  at  z = ±h/2  are  to  be  kept  at  while 
the  surface  temperature  will  be  raised  from  T\  to  T^.  Assuming  that  the  mass 
of  moisture  contained  in  the  voids  per  unit  volume  of  void  space  on  the  bound- 
ary is  constant,  then  the  relative  humidity  of  the  ambient  air  will  decrease  as 

★ 

the  temperature  is  raised  . The  opposite  occurs  when  the  RH  is  increased.  As 
it  is  to  be  expected,  an  increase  in  the  ambient  temperature  will  cause  moisture 
desorption  while  a decrease  in  the  ambient  temperature  leads  to  moisture  absorp- 
tion. These  results  are  summarized  in  graphical  form  for  the  T300/5208  epoxy 
resin  with  the  coupling  constants  u = 0.1,  X = 0.5  and  v = 0.5  as  determined 
earl ier. 

Figures  10  to  12  give  plots  of  m(t)  as  a function  of  2z/h  for  Ti  = 21 °C  and 
(RH)i  = 52"  and  three  different  values  of  Tf  * 0°C,  10°C  and  12.78°C.  Each 
graph  contains  6 curves  corresponding  to  the  time  elapsed  starting  from  t = 11.7 
hrs  to  1,961  hrs  inclusive  as  indicated.  For  these  cases,  AT  is  negative  and 
the  moisture  level  in  the  plate  will  increase  with  time.  The  differences  be- 
tween the  curves  diminish  as  the  equilibrium  condition  is  approached.  With  the 
same  initial  conditions,  Figures  13  to  15  display  the  results  for  AT  positive 

As  an  example,  if  initially  T^  = 21°C,  P.  = 0.18125  psi  and  (RH)..  = 52%,  then 

the  (RH)f  of  the  final  state  for  Tf  = 43°C  and  P^.  = 1.2748  psi  can  be  calculated 
from 

<RH>f  ■ rf  <T7>  • rasf  <Sf)  * ,5-28* 

where  Pf  is  the  saturated  vapor  pressure  at  the  temperature  Tf. 


r 


where  = 43°C,  63°C  and  83°C  are  all  greater  than  = 21 °C.  The  opposite 
trend  is  observed,  i.e.,  the  moisture  level  in  the  plate  will  now  decrease  with 
time  until  an  equilibrium  state  is  reached.  The  influence  of  AT  on  m(t)  can  be 
best  illustrated  by  plotting  m(t)  versus  /tat  2z/h  = 0.5  as  shown  in  Figure  16. 

The  curves  for  the  6 different  values  of  T.^  offer  a quantitative  assessment  of 
moisture  absorption  and  desorption  as  AT  changes  sign. 

Influence  oi  coupling.  An  attempt  is  made  in  Figure  17  to  illustrate  the 
difference  between  the  coupled  solution  obtained  from  equations  (3)  and  the  un- 
coupled equation 

1_  (d  1£)  . = o (37) 

Note  that  for  (RH)..  = 52%,  the  dotted  curves  based  on  equation  (37)  can  differ 
significantly  from  the  solid  curves  of  the  coupled  theory.  In  the  case  of  (RH)i 
= 100%  and  T^  = 21°C  changes  to  Tf  = 43°C,  the  difference  is  even  more  appreciable. 
Figure  18. 

The  foregoing  results  reveal  that  the  coupling  of  moisture  and  heat  is  in- 
herent in  the  study  of  sudden  temperature  change  in  composites  at  a given  moisture 
level.  The  extent  to  which  coupling  influences  the  mechanical  behavior  of  the 
composite  can  be  evaluated  by  calculating  for  the  stresses  and/or  strains.  This 
will  be  done  in  the  section  to  follow. 

TRANSIENT  STRESSES 

The  mechanical  behavior  of  composites  may  be  altered  when  exposed  to  high 
temperature  and/or  moisture  environments.  Their  behavior  should  be  understood 
before  the  full  potential  of  composites  can  be  realized.  For  the  resin-base  com- 
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Ipos i te  treated  earlier,  moisture  is  assumed  to  diffuse  into  the  solid  in  much 
the  same  way  as  heat.  It  tends  to  degrade  the  mechanical  properties  and  intro- 
duce dimensional  changes  of  the  composite  similar  to  those  caused  by  thermal 

i 

changes.  In  this  preliminary  analysis,  the  coupling  of  these  two  effects  will 
be  assumed  only  in  the  diffusion  process  while  the  hydroelastic  and  thermal  elas- 
tic stresses  are  taken  to  be  additive. 

Bai-tc  zquationA.  Let  the  plate  in  Figure  1 extend  to  infinity  in  both  the 

★ 

x and  y directions  and  be  free  from  mechanical  loads  . The  material  of  the  plate 
is  assumed  to  be  isotropic  and  homogeneous.  Only  hydrothermal  stresses  will  be 
treated,  i.e., 

V a.j  = E(e.j  - aAT  - SAm)  (38) 

I 


where  E is  the  Young's  modulus  of  elasticity,  a the  coefficient  of  thermal  expan- 
sion and  S the  coefficient  of  moisture  expansion.  The  quantity,  m,  is  defined 
by  equation  (19).  Since  the  stress  state  is  a function  of  the  thickness  vari- 
able z only,  shear  stresses  vanish  everywhere  and  a. j and  e^.  consist  of  normal 
components  only.  Hence,  the  strain  and  stress  relations  may  be  written  as 

s ■ r [VVP  lVa2>] 

zy  ' I fVvP  (V‘z^  (39) 

£z  ■ I tvvp  (ox%^ 

★ 

The  stresses  induced  by  mechanical  loads  can  simply  be  added  into  those  due  to 
moisture  and  heat. 
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where  v is  the  Poisson's  ratio  and  a,  = 0 will  be  assumed. 

P z 

For  the  isotropic  and  homogeneous  material,  the  stresses  induced  by  the 
strains 


e°  = = - oAT  - SAm 

x y 


are  given  by 


o o _ aEAT  _ SEAm 

ax  ~ ay  ' 1 “V  * l-v„ 
3 P P 


which  prevails  everywhere  in  the  plate. 


In  order  to  free  the  plate  edges  from  external  stresses,  it  is  necessary  to 
apply  stresses,  equal  in  magnitude  and  opposite  in  direction,  to  those  of  equa- 
tion (41).  The  following  average  tensile  stresses 


a = a = 
x y 


[aE  / ATdz  + SE  / Amdz] 


are  thus  introduced.  The  final  result  for  a free  edge  plate  is 


°x  ■ ■ T^T  <T-T>  * »-"> 

J p p 


in  which  T and  m are  the  temperature  and  moisture  averaged  through  the  plate 
thickness . 

change..  Based  on  the  diffusion  results  obtained  earlier,  equation 
(43)  is  applied  to  find  the  stresses.  For  the  T300/5208  resin,  the  following 
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material  properties  are  used: 


a = 45  x 10-6  cm/cm/°C  (25  x Id-6  cm/cm/°F) 


S = 2.68  x 10~J  cm/cm/%  m(t)  H2O 


vp  - 0.34 


E = 3.45  Gpa  (0.5  x 106  psi)  at  21°C 


The  hydrothermal  stress  distribution  throughout  the  plate  thickness  is  shown 
graphically  from  Figures  19  to  23  as  RH  is  changed  from  0%  to  13%,  33%,  52%,  75% 
and  98%  with  T..  = = 21°C.  Initially,  both  moisture  and  temperature  are  at 

the  equilibrium  state  and  hence  give  rise  to  no  stress.  As  the  relative  humidity 
on  the  plate  surfaces  is  altered,  moisture  absorption  begins.  This  causes  con- 
traction and/or  expansion  of  the  material  elements  and  leads  to  hydrothermal 
stresses  that  vary  as  a function  of  z and  t.  It  can  be  easily  seen  from  the 
graphs  that  the  stresses  near  the  surface  are  compressive  and  their  magnitude 
decrease  as  the  plate  thickness  is  increased.  These  stresses  become  tensile  in 
regions  close  to  the  center  of  the  plate.  They  increase  in  magnitude  reaching 
the  maximum  value  at  t = 1,571  hr  and  then  begins  to  decrease  settling  at  the 
zero  equilibrium  state.  This  trend  is  similar  to  the  results  of  the  uncoupled 
theory  [2]  »nd  is  to  be  expected  since  the  influence  of  coupling  due  to  diffu- 
sion was  weak  for  the  case  of  sudden  moisture  change. 

The  variations  of  the  stresses  a (or  a ) for  z=0  and  ±h/2  with  time  are  sum- 

x y 

marized  in  Figures  24  and  25,  respectively.  Figure  24  shows  clearly  that  the 
stresses  at  the  midplane  are  tensile.  They  rise  quickly  to  a peak  and  then  de- 
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cay.  Their  amplitude  increases  with  the  relative  humidity  of  the  environment. 
The  time  variation  of  the  compressive  stresses  at  the  plate  surfaces  is  similar 
except  that  the  peaks  are  much  more  pronounced.  This  is  illustrated  in  Figure 
25. 


Tempesicrfute  change,.  Now,  let  the  moisture  concentration  on  tf  plate  be  a 
constant  and  change  the  surface  temperature  which  is  initially  kept  at  the  am- 
bient condition  T.  = 21°C  and  (RH)i  = 52".  Figures  26  to  28  show  the  results 

of  a (or  a ) against  2z/h  for  Tf  = 0°C,  10°C  and  12.78°C.  This  corresponds  to 
x y t 

a temperature  drop.  The  tensile  stresses  in  the  interior  increase  in  magnitude 
while  the  compressive  stresses  near  the  plate  surface  decrease  in  magnitude. 

The  peak  tensile  stress  occurs  at  t - 1,961  hr.  The  opposite  trend  is  observed 
when  the  surface  temperature  is  raised.  Figures  29  and  30  give  the  results  for 
T^  = 43°C  and  63°C.  The  stresses  at  the  center  region  now  becomes  compressive 
and  those  near  the  surface  are  tensile.  Maximum  value  of  the  compressive  stress 
at  z=0  occurs  at  t = 791  hr. 

The  time-dependent  character  of  the  stresses  is  exhibited  in  Figures  31  and 

32  for  z=0  and  z = ±h/2.  For  z=0,  a (or  a ) increases  in  amplitude  to  a peak 

x y 

and  then  decreases  for  negative  AT  while  a (or  a ) attains  an  oscillatory  char- 

x y 

acter  when  AT  is  positive.  On  the  surface  where  z = ±h/2,  all  the  stresses, 
whether  tensile  or  compressive,  reach  a peak  and  then  reduce  to  the  equilibrium 
condition  of  zero  stress. 

Computlion  uiith.  uncoupled  theony.  The  stress  results  for  the  uncoupled  the- 
ory are  also  observed  such  that  a comparison  with  the  coupled  theory  can  be  made 
Figures  33  to  39  for  the  uncoupled  case  correspond,  respectively,  to  the  results 
in  Figures  26  to  32  for  the  coupled  case.  Although  the  general  trend  of  the 


curves  may  be  similar,  there  are  noticeable  differences  in  the  stress  amplitudes. 
n In  order  to  be  more  specific.  Table  1 gives  a comparison  of  the  stresses  at  z=0 

and  z = ±h/2  for  6 different  values  of  the  final  temperature.  The  percent  of 
deviation  between  the  results  of  the  coupled  and  uncoupled  theory  is  calculated 
for  elapsed  time  t = 11.7  hr,  402  hr  and  1,961  hr.  The  largest  deviation  occurs 
at  t = 1,961  hr.  Note  that  for  positive  AT,  i.e.,  temperature  increase,  the 
coupling  of  moisture  and  heat  can  alter  the  stress  anywhere  from  40  to  80  percent 
depending  on  AT.  In  such  cases,  the  stresses  predicted  from  the  uncoupled  theory 
may  not  adequately  model  the  physical  problem. 

CONCLUSIONS 

For  the  T300/5208  epoxy  resin  composite  material  treated  in  this  study,  it 
is  seen  that  the  interaction  of  moisture  and  heat  can  significantly  alter  the 
stress  distribution  in  the  composite.  Although  thermal  and  moisture  diffusion  do 
not  peak  simultaneously  because  of  the  wide  margin  of  difference  between  the  co- 
efficients a and  6,  the  way  in  which  moisture  and  heat  interact  in  a solid  is 
complicated  and  cannot  be  disposed  on  intuitive  grounds.  In  particular,  when  a 
composite  is  subjected  to  a sudden  temperature  change  on  its  surface  the  tran- 
sient stresses  predicted  from  the  coupled  and  uncoupled  theory  can  differ  appre- 
ciably. 

In  this  preliminary  analysis,  material  isotropy  and  homogeneity  have  been 
assumed.  These  simplifications  should  be  further  investigated  by  incorporating 
the  real  structure  of  the  composite.  What  lies  ahead  is  the  formulation  of  a 
finite  difference  method  that  treats  three  independent  variables:  two  in  space 
(x,y)  and  one  in  time  (t).  This  will. permit  an  evaluation  on  the  effect  of  ma- 
terial anisotropy,  the  presence  of  cavities  and  nonuniform  temperature  and/or 
moisture  boundary  conditions.  These  additional  influences  will  also  interact 
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TABLE  1 - PERCENT  DEVIATION  OF  STRESSES  IN  T300/5208  FOR 
COUPLED  AND  UNCOUPLED  THEORY  WITH  T.  = 21 °C  AND 
(RH).  = 522  1 


Stress  at  Surface  Stress  at  Midplane 

(psi)  (psi) 


Final 

Temperature 

Time  (hr) 

Coupled 

Uncoupled 

Dev.  2 

Coupled 

Uncoupled 

Dev.  ' 

0°C 

11.70 

-1622.82 

-1637.52 

0.91 

85.46 

70.75 

17.2 

401.5 

-1685.47 

-1931.90 

14.63 

196.41 

137.76 

29.9 

1961.0 

-1566.77 

-1966.08 

25.49 

436.03 

536.38 

23.0 

10°C 

11.70 

-1015.07 

-1024.26 

0.91 

56.13 

46.94 

16.37 

401.5 

- 989.52 

-1134.84 

14.69 

166.18 

139.11 

16.29 

1961.0 

- 837.88 

-1060.13 

26.53 

345.44 

427.36 

23.71 

12.78°C 

11.70 

- 608.12 

- 613.64 

0.91 

31 .77 

26.25 

17.37 

401.5 

- 607.79 

- 696.35 

14.57 

107.10 

89.68 

16.26 

1961.0 

- 509.24 

- 646.47 

26.95 

229.71 

285.67 

24.36 

43°C 

11.70 

402.05 

405.88 

0.95 

- 0.96 

2.87 

401.5 

513.22 

586.51 

14.28 

-203.48 

-192.23 

5.53 

1961.0 

226 . 98 

317.31 

39.8 

-133.22 

-184.71 

38.65 

63°C 

11.70 

51 .81 

52.78 

1.87 

46.97 

47.94 

2.06 

401.5 

358.10 

420.71 

17.48 

-238.39 

-252.01 

5.71 

1961.0 

47.35 

67.25 

42.03 

- 28.15 

- 39.39 

39.92 

93°C 

11.70 

-465.37 

-468.67 

0.71 

106.22 

102.92 

3.11 

401.5 

64.35 

113.84 

76.9 

-138.60 

-190.21 

37.24 

1961.0 

- 5.28 

- 0.84 

84.1 

2.82 

0.299 

89.40 
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with  moisture  and  heat  and  should  be  assessed  quantitatively  such  that  their 
individual  contribution  on  the  overall  mechanical  behavior  of  the  composite  can 
be  understood. 
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Figure  3 - Variations  of  moisture  content  with  time  caused  by 
sudden  change  in  moisture  for  T300/5208 
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Figure  4 - Sudden  moisture  change  from  RH  = 0%  to  13%  at  21 °C 
for  T300/5208 
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Figure  5 - Sudden  moisture  change  from  RH  = 0%  to  33%  at  21 °C 
for  T300/5208 
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Figure  9 - Moisture  absorption  speed  for  T300/5208  in  RH  = 100% 

air  with  different  temperature  based  on  uncoupled  theory 
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Figure  15  - Sudden  temperature  change  from  21 °C  (RH  = 52") 
to  83°C  for  T300/5208 
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Figure  17  - Comparison  of  coupled  and  uncoupled  results  for  moisture 
absorption  and  desorption  speed  with  ambient  condition  of 
21 °C  and  RH  = 52%  (T300/5208) 
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Figure  18  - Moisture  desorption  due  to  temperature  change  from  21°C 
CRH  = 100%)  to  43°C  for  T300/5208 
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Figure  26  - Stress  variations  for  temperature  drop  with  T^.  = 0°C 
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Figure  35  - Uncoupled  stress  solution  for  AT  = -8.22°C 
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